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Analysis  of  the  Spectral 
Signature  of  Breaking  Waves 

P.A.  Hwang  and  D.W.  Wang 
Oceanography  Division 

Introduction:  Wave  breaking  plays  an  important 
role  in  air-sea  interaction  processes  that  affect  the 
worlds  climate  and  the  cycle  of  greenhouse  gases 
in  Earths  atmosphere.  Over  many  decades,  wave 
breaking  is  assumed  to  be  spectrally  broadbanded. 
Investigation  of  source  function  balance  of  short  ocean 
waves,  however,  leads  to  the  conclusion  that  dissipa¬ 
tion  function  displays  a  quasi-singularity  behavior.  In 
other  words,  the  spectral  property  of  wave  breaking 
is  localized  in  wavenumber,  with  a  strong  signature  at 
wave  components  near  1  m  or  shorter,  depending  on 
the  degree  of  swell  influence.  The  result  is  consistent 
with  observed  properties  of  radar  sea  spikes  and  fre¬ 
quency  characteristics  of  sea  surface  sound  generated 
by  bubble  clouds. 

Analysis:  For  wind-generated  waves,  the  dissipa¬ 
tion  source  term  remains  the  most  difficult  to  formu¬ 
late.  Reference  1  describes  an  approach  to  establish 
the  dissipation  function  in  the  equilibrium  region  of 
gravity  waves  where  nonlinear  wave-wave  interaction 
is  much  weaker  than  wind  input  and  breaking  dis¬ 
sipation.  Theoretical  and  experimental  studies  of  the 
energy  transfer  from  wind  to  waves  have  led  to  the 
wind  input  parameterization  function, 

SJJk)  =  mo  ( ujc )2  N{k),  (1) 

where  k  is  wavenumber,  m  **  0.04,  a  is  intrinsic  fre¬ 
quency,  is  wind  friction  velocity,  c  is  wave  phase 
speed,  and  N  is  wave  action  density.  Expressed  in 
terms  of  dimensionless  spectral  density,  B  =  k*  =  oNIg, 
where  g  is  gravitational  acceleration,  the  input  and  dis¬ 
sipation  functions  can  be  written  as  Sw{k)  =  m(ujc )2 
gk~4B(k)  and  D  =  gk^f  \B).  Equating  Sw  and  D  leads 
to  the  implicit  equation  for  the  unknown  dissipation 
relation  f{B), 
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where  f;  =  ( ujc )2.  It  is  concluded  that  if  the  variation 
of  B  with  uj c  can  be  established  reliably,  f  \B)  can  be 
derived  from  Ref.  2. 


Application  to  Field  Data:  By  applying  a  free 
drifting  operation  to  minimize  Doppler  frequency 
effect,  spectra  of  short  waves  between  0.02  and  6  m 
(0.5  to  12  Hz)  are  collected  from  the  open  ocean.2 
The  environmental  conditions  cover  a  wide  range 
of  wind  speeds  and  sea  states  (Table  1).  The  data  are 
further  divided  into  wind  seas  (Subset  1)  and  mixed 
seas  (Subset  2).  Figure  4  shows  scatter  plots  of  B(ujc); 
the  data  can  be  represented  by  a  power  law  function, 
B(ujc )  =  A0(ujc)a°  =  A0^a°/2.  The  coefficient  and 
exponent,  A0  and  a0,  are  given  in  Figs.  5(a)  and  5(b). 
For  a  power  law  dissipation  relation,1  f(B)  =  AjBad,  A ^ 
and  a ^  can  be  calculated  by 

<^=1+  —  ,  (3a) 

Ai 

Ad=may)“ad.  (3b) 

The  results  are  given  in  Figs.  5(c)  and  5(d). 

Discussions  and  Conclusions:  The  presence  of 
a  localized  region  in  wavenumber  where  the  spectral 
density  is  only  weakly  dependent  on  the  forcing  wind 
condition  (Fig.  5(b))  suggests  that  wave  growth  is 
quenched  by  strong  dissipation  in  that  region.  For 
Subset  1  (wind  waves),  a0  decreases  to  less  than  0.5, 
with  a  minimum  of  0.22,  for  X  between  0.16  and  2.1 
m.  Correspondingly,  the  dissipation  function  in  this 
region  shows  a  very  nonlinear  dependence  on  B,  as 
strong  as  Bu)  and  approaches  a  delta-function  (Fig. 
5(c)-(d)).  These  wavelengths  are  sufficiently  long  that 
viscous  dissipation  or  parasitic  capillary  generation 
as  an  energy  sink  of  gravity  waves  may  be  neglected. 
Wave  breaking  as  the  major  dissipation  mechanism 
can  be  reasonably  assumed.  These  results  indicate  that 
the  length  scale  of  breaking  wind  waves  is  localized 
in  the  range  of  about  0.16  to  2.1  m.  The  presence 
of  swell  broadens  the  breaking  scales  toward  shorter 
waves  (breaker  size  ranges  between  0.2  and  0.9  m), 
possibly  through  wave-current  interaction. 
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Table  1  —  Summary  of  Experimental  Conditions 


No.  of  Spectra 

U to  (m/s) 

«*  (m/s) 

Hs  (m) 

77  (s) 

Subset  1 

291 

3.6-  14.2 

0.11-0.55 

0.21  -2.84 

1.93-6.33 

Subset  2 

106 

2.6-10.2 

0.063  -  0.40 

0.21-0.75 

2.07-4.68 
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FIGURE  4 

Scatter  plots  of  B(u*/c)  for  different  spectral  wave  components.  The  wavenumber  is  marked  at  the  top  left 
of  each  panel.  Circles  and  solid  curves  are  for  Subset  1 ,  dots  and  dashed  curves  are  for  Subset  2  (from 
Ref.  2). 


FIGURE  5 

(a]  Aq,  and  (b)  ao  of  the  empirical  power 
law  function  B(u*/c),  and  (c)  Aj,  and  (d) 
aj  of  the  dissipation  power  law  function 
f(B).  Note  that  the  range  of  Aj  covers  1  8 
orders  of  magnitude.  For  wind  seas,  the 
coefficient  increases  about  four  orders  of 
magnitude  between  4  <  k  <  6  (breaker 
size  ranges  between  1  and  1 .5  m).  In 
the  presence  of  swell,  the  breaking  scale 
moves  to  higher  wavenumber,  7  <  k  <  30 
rad/m  (breaker  size  ranges  between  0.2 
and  0.9  m).  (From  Ref.  2.) 


1 84  2005  NRL  Review  ocean  science  and  technology 


